Quantum trajectory equation for multiple qubits in circuit QED: Generating 

entanglement by measurement 
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In this paper we derive an effective master equation and quantum trajectory equation for mul- 
tiple qubits in a single resonator and in the large resonator decay limit. We show that homodyne 
measurement of the resonator transmission is a weak measurement of the collective qubit inversion. 
As an example of this result, we focus on the case of two qubits and show how this measurement 
can be used to generate an entangled state from an initially separable state. This is realized without 
relying on an entangling Hamiltonian. We show that, for current experimental values of both the 
decoherence and measurement rates, this approach can be used to generate highly entangled states. 
This scheme takes advantage of the fact that one of the Bell states is decoherence-free under Purcell 
decay. 



I. INTRODUCTION 

In the last few years circuit quantum electrodynamics 
(circuit QED) [1, 2. 3], a solid state analog of cavity QED 
[4, 5, 6], has grown into a mature field. This system is 
based on superconducting qubits [7] acting as artificial 
atoms and a distributed [1, 2, 8, 9, 10, 11, 12, 13, 14, 
15, 16, 17] or lumped [3, 18] resonator which acts as 
a harmonic oscillator. Examples of its success are the 
observation of the particle-like nature of microwave pho- 
tons [8] , generation of a single photon [9] and Fock states 
[10, 11], observation of Berry's phase [12], use of a quan- 
tum bus to couple qubits [13, 14], single artificial-atom 
lasing [15], and observation of the fundamental limit that 
exists between measurement and dephasing [16, 17]. 

Evolution of this system obeys a master equation 
(ME) which has an interaction described by the Jaynes- 
Cummings Hamiltonian [19] and decoherence processes 
which are well described by Markovian environments [20] . 
Measurement in this system is done by operating in the 
dispersive limit (where the detuning between the res- 
onator and the qubit is much larger then their coupling 
strength). In this limit, the interaction induces a qubit- 
state dependent frequency shift on the resonator. By 
measuring the resonator output voltage with a homodyne 
measurement, information about the qubit state is ob- 
tained. This measurement is a weak continuous-in-time 
measurement of the quadrature of the resonator and the 
evolution of the conditional state is described by a quan- 
tum trajectory equation (QTE) [20, 21]. The presence 
of the resonator can make the Hilbert space needed for 
simulation of the ME or QTE quite large and impractical 
(especially for a many qubit system). 

It has been shown by us in previous work [21] that 
for a single qubit a polaron transformation can be used 
to eliminate the resonator dynamics from both the ME 



and QTE. The resultant ME and QTE have an extra de- 
cay channel which represents measurement induced de- 
phasing. Moreover, measurement is found to be a weak 
measurement of the qubit inversion operator cr^. This 
transformation has the advantage of being exact for the 
average evolution. Extending to more than one qubit 
is non-trivial and will be discussed in a future publi- 
cation [22]. Here we follow a simpler approach which 
yields very accurate results in the large resonator damp- 
ing case. This large damping limit is particularly useful 
because it corresponds to a good qubit measurement. We 
use an adiabatic approximation, similar to that presented 
in Refs. [23, 24, 25, 26], to obtain an effective ME and 
QTE. From these equations, we find that in this limit 
homodyne measurement of the resonator corresponds to 
a measurement of the collective qubit inversion (^^ ^i^l , 
where 5i is a dimensionless parameter determined by the 
system parameters). 

As an application of the derived many-qubit QTE, we 
consider the case of two qubits and tune the system pa- 
rameters such that 5i = 62- In this situation, a mea- 
surement has three possible outcomes. For one of these 
outcomes, measurement cannot distinguish between an 
excitation being in either qubit. As a result, starting with 
a particular separable input state this measurement will 
generate a maximally entangled state. Similar ideas were 
used with trapped ions to generate entanglement [27, 28] . 
Theoretical work towards realization of these ideas in cir- 
cuit QED was already done in Ref. [29]. Here however, 
we consider realistic decoherence rates as measured in 
recent experiments [30]. Moreover, to obtain large con- 
currences our scheme takes advantage of the fact that 
one of the Bell states is decoherence-free under Purcell 
decay [31]. 

The paper is organized as follows. In Sec. II we derive 
the underlying ME for many qubits coupled to a common 
resonator and adiabatically eliminate the resonator de- 
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grccs of freedom. In Sec. Ill we derive the corresponding 
QTE. In Sec. IV we investigate how entanglement can be 
generated by measurement and show that, with reaUstic 
parameters, high concurrence can easily be reached. We 
summarize our findings in Sec. V. 



II. THE MASTER EQUATION 




Hamiltonian Hctf = U'''_ffU is, to second order in A,- 



3>k 



E 



(2.4) 



where Xj = 9j/^j parametrizes the strength of the ac- 
Stark shift (fourth term) and Lamb shift (fifth term) on 
the j"^ qubit transition frequency. Jjk represents the 
strength of the coupling between the j*'' and k^^ qubit 
by virtual photons in the resonator and is given by 



J 



gj-gfc(Aj- + Afc) 



(2.5) 



FIG. 1: A schematic illustration of multiple-qubit circuit 
QED. Superconducting qubits (gray; six qubits are illus- 
trated) are fabricated inside a transmission-line resonator 
(blue). The fuU-wavelenght mode of the resonator is illus- 
trated (pink). 

We consider multiple superconducting qubits coupled 
to a transmission line resonator acting as a simple har- 
monic oscillator. This system is illustrated schematically 
in Fig. 1. In the limit where direct capacitive qubit- 
qubit coupling can be ignored, the system is described by 
the niulti-qubit Jaynes-Cummings Hamiltonian [19, 32] 
(n = l). 



(2.1) 



In this expression, ojr is the resonator frequency, ojq^ 
the j*''-qubit transition frequency, and gj the j*'^-qubit- 
resonator coupling strength. 8 ^ 8m + £c represents ex- 
ternal driving of the resonator, with £„i referring to the 
measurement drive (at frequencies close to the resonator) 
and £c referring to the control drives (at frequencies close 
to the qubits). We take 



(2.2) 



with uj,„ being close to uJr and where, for all fc, uji^ is far 
detuned from the resonator. 



In the dispersive regime, where |Aj 



> 



\gj\, we can adiabatically eliminate the resonant Jaynes- 
Cummings interaction using the transformation [41] 



U : 



exp 



[^A,(at.- 

j 



(2.3) 



where Xj ~ gj/^j is a small parameter. The effective 



Coupling to additional uncontrolled degrees of freedom 
leads to energy relaxation and dcphasing in the system. 
Integrating out these degrees of freedom leaves the qubit 
plus resonator system in a mixed state p{t) whose evolu- 
tion can be described by the ME 



p = - i[i?cfr, p] + KV[a]p + ^ li,]'D[aj ]p 

3 



(2.6) 



where 



V[A]p^ApA^-{A^A,p]/2. 



(2.7) 



In this expression, k, is the resonator decay rate, 71 j , 7<^j 
represents relaxation and dephasing on the qubit and 
the last term represents correlated relaxation due to the 
Purcell efi'ect [31, 33, 34]. While this terms is of order 
A^, it is kept as in the following adiabatic approximation 
we will require k to be large such that the product A^k 
is not necessarily small. 

To derive an effective ME for the qubits only, we start 
by removing the fast dynamics of the resonator. This 
is done by moving to the interaction frame rotating at 
the measurement drive frequency, u)m and by making the 
standard rotating- wave approximation. This allows us to 
rewrite Eq. (2.4) as 

iJcff = A^,ata + ^Xja^aCTj + (e;,a + e™a^) + i7,, (2.8) 



where A,. = lj,. — lj„i and Hq represents the Hamiltonian 
of the isolated qubits 



a|+^J,fe(a^. a++a+a,) 



E 

3 



3>k 

a7). 



(2.9) 
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Next, we move to the frame defined by 



p°(i)-DtHp(t)D[a] 



(2.10) 



where D[a] = exp[aa^ — a*a] is the displacement opera- 
tor. Applying this to Eq. (2.6) and choosing 



(2.11) 



k/2 



yields 



j 

(2.12) 



where C is the Lindblad superoperator representing only 
qubit dynamics and is given by 



£p = - i[H^,p] - ix\a\^ Y.^S,a],p] + ^7i.,25K']/ 



j:l^V[a^]p + .V 



(2.13) 



In this equation the second term represents the ac- 
Stark shift on the qubit transition frequency, with x ~ 
J2j Xj/^ and N being the total number of qubits. Sj = 
Xj/x is the fractional effect that Xi gives to the average. 

Following Refs. [23, 24], we now make an adiabatic 
approximation. That is, we will assume that quantum 
fluctuations in the displaced resonator state are small. 
In this situation, it is reasonable to assume that matrix 
elements p„„n with n, m being the resonator photon num- 
ber, decay rapidly with increasing ri, m. As a result, we 
will assume that there is a small parameter s such that 
Pnm oc e"+™ [23, 24]. Expanding the total state matrix 
to second order in e we find 

P° =Poo|o)(o| + pio|i)(o| + m|o)(i| + pii|i)(i| 

+ P2o|2)(0|+po2|0)(2| + 0(£3) 



such that the reduced state for the qubits is given by 
Q = Tr[p°] = Poo + pn. 

Substituting this expansion in the above ME leads to 



the following coupled differential equations 

Poo =-^^00 + npn + i^XjiapoiCTj - a*cr|pio), 

j 

pio =Cpio - Kpio/'2 + i^axj{pii<j] ~(j]poo) 

j 

~ ^^XjicTjPio + a*V2ajP2o) - iArPw, 
j 

pii - Kpii +i'^X]{o!*piocr] - ceajpoi) (2-15) 

j 

j 

ha =^P2Q - np2o - i2ArP20 

- i Xj (av^cr|Pio + 2cr|p2o)- 



By looking closely at these expressions, we find that the 
higher order terms are only populated at rate |a| Xj/'* 
and decay at rate k. Thus, taking 



(2.16) 



as our small parameter, we can assume that the off- 
diagonal terms pio and p2o decay much faster than the 
diagonal terms and as such can be approximated by their 
steady-state value 



Pio 



P2a = J2 

j 



iAr + k/2 
-iaV^XjCr^jPw 

K + i2Ar 



(2.17) 



Substituting these expressions into the diagonal compo- 
nents leads to the effective ME [42] 



g~tK[{Y^S,a]f,g]. (2.18) 



In this expression, Ed is the measurement-induced de- 
phasing rate and is a resonator-induced Ising-like cou- 
pling. These are 



Ed = -ir^ 5-77 and K 

A2 + k2/4 



-2AAa\^X^ 



(2.19) 



In the limit of a single qubit, the measurement induccd- 
dephasing rate obtained here correctly agrees with the 
large k limit of Ref. [17] (In this paper, this expression is 
labeled Em) and with the large k and steady-state limit 
of the rate presented in Ref. [21]. We note that one could 
use a multi-qubit polaron transformation to get results 
valid in the small k limit and which would take into ac- 
count initial transients in the resonator [22] . 
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III. QUANTUM TRAJECTORY EQUATION 



allows us to write the effective QTE in the form 



Although direct detection of the transmitted mi- 
crowave photons is possible [9] , here we will consider ho- 
modyne processing [20] . That is, we will assume that the 
signal coming from the output port of the resonator is 
mixed with a strong local oscillator tuned to the signal 
frequency and of phase (f). Given the homodyne mea- 
surement result J{t), we can assign to the qubit and the 
resonator the conditional state pj{t) whose evolution is 
governed by the QTE [35] 

pj ^ Cpj + i^[Q^,pj]^it) + ^M[2I^]pjm, (3.1) 

with C given above. Here, M [c] is the measurement su- 
peroperator defined as 



M[c]p^{c-{c)t)p/2 + p{c-{c)t)/2, 



(3.2) 



-dependent field com- 
and 2Q0 = — z 



where (c)t = TT[cpj{t)] and the 
ponents are 2/^ = ae^"^ + a^e"^ 

la^e**^. 77 is the efficiency at which the photons com- 
ing out of the resonator are detected. For the cur- 
rent circuit QED experiments, this can be written as 
Vdct = l/(-^th + 1) with iVth being the number of noise 
photons added in the amplifier stage. The measurement 
record is 



Q.i ~iK[C^5ja^^f,Qj] 



m 

(3.8) 



with 



(3.9) 



and Fci + Fba = ^i^m- This last equality is also found in 
Ref. [21] and reflects the fact that maximum information 
about the qubit is obtained by setting the phase of the 
local oscillator such that Fba is zero [Y = 0). At this 
point, the rate of information gain is Fci = V^m a-nd for 
f] = \ this system reaches the quantum limit (Fd/Fd = 
1/2). 



IV. ENTANGLEMENT BY MEASUREMENT 



(3.3) 



where ^(t) is Gaussian white noise and represents the 
photon shot noise. It is formally defined by E[^(i)] = 
and E[^(t)^(t')] = 5{t — t'), with E denoting an ensemble 
average over realizations of the noise ^{t). 

To obtain an effective QTE for the qubits only, we 
apply the small e expansion to the above QTE. For the 
stochastic part only, we find 



Poo ^\/Krj {pwe + poie"^ - {2I^)poo) ^{t), 
Pii =V7^(-(2/0)pii)^(t), 



(3.4) 



which gives the following stochastic term to the qubit 
equation 

^ (pioe-^^ + poie"^ - {2I^)g) ^(0- (3-5) 



Using the steady state value for pio given in Eq. (2.17) 
this yields 
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(3.6) 

where X ^ 3fi[a*e*('^-*)] and Y = 3[a*e'('^-^)] with 
tan(0) = K/2Ar. Defining the measurement rate Fci, the 
extra non-Heisenberg backaction Fba and the maximum 
measurement rate Fm as 



Fci 



A2 + K2/4' 



4:m]x^Y^ 

A2+K2/4' 



and Fn 



A2 + K2/4' 

(3.7) 



1 2 3 4 5 

Time, t [Fci-i] 

FIG. 2: Three typical trajectories of the conditional state 
elements pee (blue solid line) pgg (red dotted line) and p++ 
(green dashed line) when the initial condition is the separable 
'state given by Eq. (4.2). The system parameters are Fd = 
rci/2, and 71,1 — 71,2 = 7p = 0. In panel a), the collapse of 
the wavefunction is to the pure state jee), in panel b) to \gg} 
and in panel c) to the entangled state j-l-). 

As an illustration of application of the above QTE, 
we consider the probabilistic generation of entanglement 
from a separable state without using a two-qubit unitary. 
This approach can be particularly useful when such a 
two-qubit unitary is not present or hard to implement. 
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We consider the case where only two qubits are present 
and take = 0, Si = 62 = 1- The last equality is such 
that both states \eg) and \ge) pull the resonator to the 
same frequency and are thus indistinguishable. More- 
over, the local oscillator phase is assume to be chosen 
such that Fba ~ and we take Ai = — A2. The sign 
change in A is easily realized by putting the two qubits 
at both ends of the resonator and working with an odd 
mode. Finally, to rule out any possible entanglement due 
to unitary qubit-qubit coupling the indirect interaction 
J12 between the qubits is ignored in the numerics [43] . 

With a transmon-type qubit [36], we can safely take 
the qubit dephasing rate 7^ to be small [30]. In this 
situation, Eq. (3.8) in the rotating frame becomes 

gj ^^l.lT^[(Ti]gj +"fl.2'D[(72]Qj +7pX>[crf - (72]gj 

+ \oi + ai\ Qj + vr;;x[aj + ai]Qjm, 

(4.1) 

where J = \/T^i{<Jl + erf) + ^(0 ^'^d 7p = nX^ is the 
Purcell decay rate. 

For the system in the initial separable state 

1^) = + \me) + \9)) = \{\ee) + \gg) + ^2]+)), 

(4.2) 

where 

l±) = i=(|eg)±lffe)), (4.3) 

the QTE (for 71 1 = 71 2 = 7p = 0) leads to a collapse to 
the state ]ee) with probability 1/4, \gg) with probability 
1/4 and to the entangled state ]+) with probability 1/2. 
This is illustrated in Fig. 2 where typical quantum tra- 
jectories for the elements (blue solid) pgg (red dotted) 
and P+-I- (green dashed) are plotted as a function of time 
for these three possible outcomes. 

To connect this to the language of POVM's [37], or to 
a quantity measured in an experiment, we defined the 
integrated current [21, 38] 

s{t) = f J{t')dt'. (4.4) 
Jo 

If we break s into three regimes defined by s < — fth, 
—vth < s < vth, and s > vth, where we refer to i^th as the 
threshold then the POVM elements Eg, Eq, and E^ are 
measured, respectively. These are 

Eg = Ggllg + t gllo + CgHe, 

Eo = aollg + bono + con^, (4.5) 

Ee = ttellg + bello + Ce^i-e- 

where Hg = \gg){gg\, He = \ee){ee\ and Uq = \eg){eg\ + 
\ge){ge\ are projectors and the rest of the parameters arc 
simply positive real numbers (probabilities) which satisfy 
Qg + ao + Qe = 1 , bg + bo + be = 1 , and Cg + co + Ce = 1 . For 




FIG. 3: Histogram of 10000 trajectories at time t — IF^ 
a) and time t = lOF^^ b). In both panels the full blue line 
shows the expected distribution. The center peak in panel b) 
corresponds to the entangled state j-l-) while the side peaks 
to the separable states \gg) and jee). Average concurrence 
(blue solid line) and probability of success (red dotted line) 
are plotted as a function of Vtii for times t = IF^^ c) and 
t — lOF^^ d). Other parameters are the same as in Fig. 2. 



example in the POVM element Eg, Ug is the probability 
that the measurement was of the projector we desired 
and bg and Cg arc the probability of the false positive 
events IIq and lie respectively. 

A histogram of s for 10000 trajectories is plotted in 
Fig. 3 for 71.1 — 71.2 = 7p = 0. Panel a) is taken at the 
integration time t = IF^JT^ and panel b) at t = 10F~j^. 
Clearly, the measurement at the earlier time is not pro- 
jective (the false positive rates are high), whereas at the 
later time the measurement is projective since the distri- 
butions are well separated. At that time, it is possible 
to create and distinguish the entangled state ]-|-) from 
\gg) and ]ee). This is shown more explicitly in panels 
c) and d) where the average concurrence [39] and the 
probability of success Pg (defined as the number of de- 
tections in the range [— I'th, fth] divided by total number 
of measurements) as a function of the threshold i^th are 
plotted. Panel c) corresponds to t = IF^^ and panel d) 
to t ^ 10^-^. We see that, at the shorter time, it is im- 
possible to generate the entangled state while at the later 
time a highly entangled state can be created by setting a 
low threshold. If this threshold is too small the probabil- 
ity of generating the entangled state goes to zero and if it 
is too large then the unentangled results are included in 
the post selection states. However, if the measurement 
time is long enough there can be a large range of values 
for i/th [8-12 in Fig. 3 d)] where it is possible to create 
the desired entangled state with probability 1/2. 

We now take into account realistic values for the qubits 
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FIG. 4: The average concurrence given that s is be- 
tween ±y/Vcit as a function of time for 5000 tra- 
jectories and (7p,7i,i,7i,2,rd) = (0, 0, 0, l/2)rci 
(red dotted line), (l/2,0,0,l/2)rci (green dashed 
hne), (1/2, 1/200, 1/200, l/2)rci (blue solid line), and 
(1/2, 1/200, l/200,10)rci (pink dash-dotted line). The inset 
is a histogram of 5000 trajectories at time t — 5r_^^ and 
(7p, 71,1, 71,2, Td) = (l/2,0,0,l/2)rci. 

relaxation rate. We choose 7p = rci/2 and 71.1 = 71^2 = 
rci/200 as measured in Ref. [34]. These rates correspond 
to Tp = l/7j, = 40 ns and Ti = 4 //s, and to a signal-to- 
noise ratio rci/7p of 2 consistent with experimental obser- 
vations. Numerical results using these values presented in 
Fig. 4 show that states that are close to maximally entan- 
gled can be obtained even in the presence of qubit decay. 
This is because the Purcell effect, which is the dominat- 
ing source of decay for the transmon [34], docs not act 
on the entangled state [-I-), since (erf — o'2^)|+) =0. This 
state is thus a decoherence-free subspacc with respect to 
this decay channel. It is also worth pointing that it is an 
eigenstate of the Hamiltonian, such that the prepared [+) 
states are robust against further evolution of the system. 

To maximize the amount of entanglement generated, 
we find that it is best to use an integration time which is 
larger then l/7p (so that errors due to to [ee) decaying 
into [— ) have subsided) but shorter than the single qubit 
relaxation time 1/71,1(2) (causing errors with [+) decay- 



ing into \gg))- Doing this, we find that the maximum 
attainable concurrence can be as large as 0.94. Inter- 
estingly, these results depend on the signal-to-noise ratio 
rci/7p, but not on the efficiency of the detector 77 which 
was so far taken to be unity. To show this, we have 
simulated the average concurrence for s between ztyT^t 
as a function of time and for an efficiency 77 = 1/20 
(Fd = lOFci). This is shown in Fig. 4c) as the pink 
dash-dotted line. These results are equal, within statis- 
tical error, to the situation where 77 = 1 (blue solid line). 
This is simply because the extra dephasing caused by 
the inefficiency only results in the system losing coher- 
ences between the measurement eigenstates faster than 
the limit imposed by measurement (the quantum limit). 
This does not change the measurement statistics as, with 
dephasing, measurement selects the final state out of a 
classical mixture rather than from a quantum superposi- 
tion. 



V. CONCLUSION 

In this paper we have derived a master and quan- 
tum trajectory equation for multiple qubits in a res- 
onator. These equations are valid in the limit where 
e = Xj\o^^ I ^ ^ 1- In this limit, we find that mea- 
surement of the resonator transmission leads to a weak 
measurement of the qubit observable ^i^f . As an ex- 
ample of this result, we have discussed how entanglement 
can probabilistically be generated by measurement only. 
We have shown that with current experimental decoher- 
ence rates, this process can yield a highly entangled state. 
This is mostly due to the fact that, in circuit QED, a ma- 
jor source of relaxation (the Purcell effect) does not affect 
one of the Bell states which is therefore protected from 
decay. 
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